X-ray method to study temperature-dependent stripe domains in MnAs/GaAs(001) 
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MnAs films grown on GaAs (001) exhibit a progressive transition between hexagonal (ferromagnetic) and 
orthorhombic (paramagnetic) phases at wide temperature range instead of abrupt transition during the first-order 
phase transition. The coexistence of two phases is favored by the anisotropic strain arising from the constraint 
on the MnAs films imposed by the substrate. This phase coexistence occurs in ordered arrangement alternating 
periodic terrace steps. We present here a method to study the surface morphology throughout this transition by 
means of specular and diffuse scattering of soft x-rays, tuning the photon energy at the Mn 2p resonance. The 
results show the long-range arrangement of the periodic stripe-like structure during the phase coexistence and 
its period remains constant, in agreement with previous results using other techniques. 



The integration of magnetic films in semiconductor de- 
vices represents one of the major challenges in materials sci- 
ence. MnAs is a very promising material for spin injec- 
tion, in combination with III-V semiconductor compounds 
(such as GaAs) 1,2 . Around room temperature (RT), bulk 
MnAs undergoes a first order structural and magnetic phase 
transitions^ between the low temperature a phase, hexag- 
onal and ferromagnetic, and the high temperature P phase, 
orthorhombic and paramagnetic. MnAs films grown on 
GaAs(OOl) 6 ' 7 ' 8 exhibit a phase coexistence over a tempera- 
ture range that can extend between °C and 55 °C, depend- 
ing on the thickness of the film^liLmi. Several local probe 
techniques have been used to characterize the surface mor- 
phology and magnetic structure, such as atomic force mi- 
croscopy (AFM) 12 , magnetic force microscopy (MFM)i^and 
X-ray magnetic circular dichroism photo-emission electron 
microscopy (XMCDPEEM)^. They highlighted the forma- 
tion of periodic terraces and magnetic domains at intermediate 
temperatures for the phase transition. Information on the long 
range order of terraces is usually gained through scattering ex- 
periments. In the case of MnAs, x-ray diffraction experiments 
showed the coexistence of the two crystallographic phases^ii, 
but could not address the formation of periodic domains with 
long range order. 

In this work we have used specular and diffuse x-ray reflec- 
tivity to follow the formation and evolution with temperature 
of the terrace morphology throughout the a — p phase transi- 
tion. In contrast to the microscopy technique, such as atomic 
force microscopy, which probes small area and only near the 
surface, the x-ray reflectivity technique proposed here ana- 
lyzes over a much larger area of the surface and deepness, 
providing long range order information and also the interface. 
On MnAs films analyzed using this technique we observe the 
formation and disappearance of diffuse satellite peaks associ- 
ated with the terrace periodicity as a function of temperature. 
The modeling of their intensity profiles yields numbers for the 
temperature-dependent terrace widths. 

A 130 nm thick MnAs was grown by molecular beam epi- 
taxy on a GaAs(OOl) substrate at 250 °C in a growth condition 



to obtain A-type orientation , i.e., the MnAs plane in the a 
phase parallel to the GaAs (001) plane with the MnAs [0001] 
c-axis along the GaAs direction. 

We have first used AFM to identify the terrace-like steps 
discussed before. Figure la depicts a room temperature AFM 
image of the surface topography of the MnAs film, which 
shows the formation of stripe-like domains, elongated along 
the (0001) direction, during the coexistence of two phases and 
a clear alternate periodic structure between a and p phases is 
observed. Over the extension of the AFM image, these ter- 
race steps appear periodic, with a modulation period of about 
600 nm. The lower terraces correspond to regions where no 
magnetic signal is obtained with Magnetic Force Microscopy 
(MFM), therefore they can be associated to the high temper- 
ature paramagnetic p phase. The higher terraces (ferromag- 
netic a phase) exhibit a complex magnetic structure where 
domains of opposite orientation seem to be intercalated^^. 
Several authors have already addressed the magnetic domain 
structure of these films using MFM techniques, and at least 
three different types of domain formations were observedi^ii. 
Here we will concentrate on our X-ray method to study the 
terrace configuration. 

In order to quantitatively evaluate the surface morphol- 
ogy of the sample, we performed resonant x-ray scattering 
measurements at beamline SU-7 of the SuperACO storage 
ring (LURE laboratory, Orsay). The beamline, equipped 
with a linear undulator source, covers the 100-1000 eV 
range with a resolving power of about 2000. The endsta- 
tion is a 2-circle (co/29) reflectometer working in ultra-high 
vacuumi^. As shown in Fig. lb, the scattering geometry was 
co-planar, with the incoming beam (of wave vector k,) im- 
pinging on 0.5mm x 0.5mm of the MnAs film at a grazing 
angle , and the scattered photons (k/) collected at an angle 
29 with respect to the incident beam. The scattering vec- 
tor q = k/ — k,- can be separated into two components, one 
parallel q x = 3? [cos (20 - co) — cos (go)] and one perpendicu- 
lar q z — y [sin (20 — co) + sin (©)] to the sample surface. A 
structure that is periodic with a modulation period I will give 
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Pi gure 1 ! a) Atomic force microscopy image of the 130 nm thick MnAs film grown on 
GaAs(OOl); b) Scattering geometry used to examine the terrace step structure. Scattering 
vectors are described in the text. 



rise to constructive interference in the scattering process when 

n - 2k 

The horizontal terrace step structure can be studied in de- 
tail by performing rocking scans of the sample, i.e. CO scans at 
fixed 20. Over a limited range in CO around specular (CO = 2 j), 
a rocking scan corresponds to a q x scan at fixed q z . In Fig. 2a 
we show the map of the scattered intensity that results from a 
series of 66 q x scans (201 points each) performed as a function 
of temperature at a given q z value (1.68 nm-1, correspond- 
ing to 20 = 30° and photon energy (640 eV). Fig. 2b shows 
the line plots of three of these scans, at temperatures corre- 
sponding to the a — p phase coexistence region. In the inter- 
mediate temperature region, we observe, besides the specular 
peak (zero order) at q x = 0, two other peaks (first order) at 
q x = ±0.0107 nm -1 . We ascribe them to the long range lat- 
eral order produced by the stripe domains in the phase coex- 
istence temperature region. The associated modulation period 
is t = — = 587 nm. It is worth noticing that, while the in- 
tensity of the first order peaks varies drastically as a function 
of temperature, their position remains unaltered. This implies 
that the modulation period of the stripes, i.e. the sum of the 
widths of the two domains with different structures, remains 
constant. The a to p transition takes place versus temperature 
with the widths of the p stripe increasing at the expense of the 
width of the a phase, their sum remaining constant. This has 
been previously observed by AFMi£ and by XMCDPEEMii. 

In order to reproduce the scattering profiles, a model was 
introduced. Following Holy et al4&, the scattering intensity 
stemming from a lateral periodic structure is given by 



Hlx^z) = const. 1 7/7/ 1 \C(q x 
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I \F(q x ,q z )\ 



where 7J and Tf are the incident and exit transmission func- 
tions, which depend on the incident and exit angles and can 
give rise to multiple scattering for very thin films or multi- 
layers (not the present film). In this case one might need to 
address the behavior of T and Tf with angle to understand the 
scattering data. 

F(q x ,q z ) EE / e- iq * h W- iq * x dx 
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Figure 2: a) Two-dimensional plot of diffuse scattered intensity measured at 26 — 30° 
as a function of q x (i.e. sample angle G)) and of temperature in the phase coexistence 
region and b) Corresponding fitted profiles following the model given in the text; c) Se- 
lected X-ray scans for different temperatures; d) Fitted profiles for corresponding scans 
on the left. Each profile was fitted separately according to equationpl 



h (x) is given by 



h(x) 



, if < x < s 
h ,ifs<x<£ ' 



where s is the width of the p phase terrace and I the sum of 
the widths of the two terraces. 

c (<?.<-) = exp \-i (R m - R n ) q x ] ) 



is the correlation function of different sets of two terraces 
averaged over the whole sample surface. This correlation 
function is proportional to a periodic sequence of 8-like 
functions 16 centered at reciprocal lattice positions at intervals 
g= 2 f> resulting in 

C(q x )=N(g £ §(q x -pg)}, 

p= -oo 

where N is the number of periods. These 8-like functions 
take a Lorentzian-like lineshape with full-width a. This finite 
width here represents the combination of the limited coher- 
ence length of the x-ray beam and the correlation length of 
the stripes. We obtain the intensity profile, which in our was 
limited to the first 4 reciprocal lattice points 



I{qx,q z ) <* \TjT f 



2 f s 2 + {e-s) 2 +2cos(q z h)£(£-s) 



q 2 x +C 2 
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1 ' P=-4(P^0) - f ) 
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is the Fourier transform of the height profile function h (x) of 
the surface, covering one terrace of each kind. In our case 



In figures 2c and 2d the results of calculations are compared 
to the corresponding experimental results of 2a and 2b. The 
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Figure 3 '. Temperature dependence of the p 1 phase terrace width. The error bars are 
estimated from the statistical fluctuation of the scattered intensities of the specular and 
first diffuse satellite peaks. The error bars are larger in the intermediate temperature 
region due to the lower intensity of the specular peak. 



Figure 4: a) Scattered intensity as a function of vertical momentum transfer q z for 
a parallel momentum transfer q z for a parallel momentum transfer q x — 0.0107nm . 
Hollow circles: T = 5°C; filled circles: T = 32°C; b) Intensity ratio /[r=32ocj+/(r=5ocj °^ 
the two curves in a), used to determine the terrace height. 



intensity profiles were fitted separately for each temperature, 
using s, I and a and an overall scaling factor as fitting pa- 
rameters. This model clearly reproduces extremely well the 
scattering data. The parameter o = (3.2 ±0.1) x \0~ A nrrT l 
remained constant throughout all measurements, which indi- 
cates that the major source of broadening comes from the fi- 
nite coherence length of the x-ray beam. 

In order to evaluate the relative weight of the a and P 
phases, one can reduce the number of fitting parameters by 
limiting the analysis to the ratio between the intensities of the 
first order peak and of the specular reflectivity. In this case the 
intensity ratio of these two peaks as predicted by equation^is 
approximately given by 17 

l(g x =f) _ 16sinfc)sin(^) 2 /(^) 2 
/ (q x = 0) s 2 + (£ - sf + 2 cos (q z h) £(£- s) 

The terrace width was determined by solving numerically 
equation independently from other fitting parameters. As 
a result, much smaller error bars for the terrace width were 
obtained. Fig. 3 shows the temperature dependence of the 
p-phase terrace width. 

Finally, we have measured the scattered intensity as a func- 
tion of q z for a constant non-zero value of q x . This can be 
obtained by independently positioning sample and detector 
angles according to a predetermined set of values. Fig. 4a 
shows the results of two q z scans (T = 5°C and T — 32°C) 
taken at q x = 0.0107 nm~ l , i.e. at q x value matching the lat- 
eral order parameter of the stripes. Besides the peak of the 
transmission function (at q z = OAnm^ 1 ) present in scans per- 
pendicular to the surface, at T = 32°C there is a second peak 
centered around q z = 1 nrrT 1 , which comes from the interfer- 
ence between the waves scattered at the lower and the higher 
terraces. 



In order to eliminate the contribution of the transmission 
function and to enhance the interference term, Fig. 4b shows 
the intensity ratio of the two scans. This ratio is centered at 
q z = 1 .08 ± 0.02 nm" 1 , corresponding to a difference in height 
between the two stripes of 5.82 nm (this height corresponds to 
an average value over the whole sample surface). 

In contrast to the surface morphologic microscopy tech- 
niques, the X-ray technique used here probes wide area and 
the depth of the film including also the different X-ray ab- 
sorption edges to select the structures involving particular 
atomic elements. This will allow the determination of chemi- 
cal/magnetic ordering of specific atoms belonging to the crys- 
tal lattice. The introduction of this technique will therefore be 
useful in special cases where microscopy techniques are not 
appropriate. 

It is also important to understand why such diffuse 
satellites were not previously observed in x-ray diffraction 
experiment*^: there, one mainly observed Bragg peaks, 
which are sensitive to atomic distances of individual terraces. 
The radiation is scattered by atoms and variations in the ter- 
race widths lead to an incoherent sum of the x-ray waves from 
neighboring terraces. In our present reflectivity method, one 
is not sensitive to the details of the crystal structure. The x-ray 
wave does not lose coherency when scattered by each terrace 
and the resulting interference pattern reveals the modulation 
of the surface structure. 

In summary, we demonstrate in this work that the resonant 
soft X-ray diffuse reflectivity method is a powerful tool to in- 
vestigate microscopic structures such as the terrace-like struc- 
tures of MnAs/GaAs(001) films. We observed a clear long 
range periodic arrangement of stripe-like micro-structures al- 
ternating between a and p phases during the phase coexis- 
tence, which takes place between 10 and 45 °C. The period of 
the modulation structure obtained is about 600 nm, which re- 
mains constant with the temperature. The results are in good 
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In this procedure we have subtracted the contribution of the spec- 
ular peak from the intensity of the first satelite peak, i.e., we have 
substracted the nearly temperature independent component of the 

i(<lx=-) 

intensity of the first satellite. In this case, ^* _y is aproxi- 
matelly given by the ratio of the second (p = 1) and first terms 
of eq. 6, leading to eq. 7. 



